Integrability, stability, and adiabaticity in nonlinear stimulated Raman adiabatic 
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We study dynamics of a two-color photoassociation of atoms into diatomic molecules via nonlinear 
Stimulated Raman adiabatic passage (STIRAP) process. This system has a famous counterpart in 
(linear) quantum mechanics, and been discussed recently in the context of generalizing quantum 
adiabatic theorem to nonlinear systems. Here we use another approach to study adiabaticity and 
stability in the system: we apply methods of classical Hamiltonian dynamics. We found nonlinear 
dynamical instabilities, cases of complete integrability, and improved conditions of adiabaticity. 



Adiabatic theorem [l[ of quantum mechanics has found 
wide applications in quantum state manipulations. Dy- 
namics of Bose-Einstein condensates (BEC) Q intro- 
duces paradigm of nonlinearity into a quantum world 
(note studies on nonlinear Landau-Zener models [H, 0], 
solitons, shock waves, etc). Along with nonlinearity, it 
also brings a question of how to analyze adiabaticity in 
nonlinear systems. It is not justified to apply the adia- 
batic condition of quantum mechanics to nonlinear BEC 
systems @Jp|. Instead, a general method was recently 
suggested Q and applied to a specific example of two- 
color Raman photoassociation system Q. Here, we con- 
sider the same example as in [y, 0] and notice that, as 
a typical nonlinear system, it possesses many counterin- 
tuitive difficulties in the analysis. The goal of our paper 
is twofold. Firstly, we (inspired by ideas of @) suggest 
another approach to analyze adiabaticity and stability in 
nonlinear BEC-related systems: we use methods of clas- 
sical Hamiltonian dynamics [Tlj . Secondly, we improve 
a theory of nonlinear STIRAP (Stimulated Raman adia- 
batic passage). 

Following [y, 0] , consider a three-level A-system shown 
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FIG. 1: (a) A three-level A-system (b) Counterintuitive STI- 
RAP sequence of pulses fl Pt d (c) Phase portrait of ([3|. As 
parameters are changing in accordance with a pulse sequence 
shown in (b) , the stable fixed point (corresponding to the dark 
state) moves from z = 0. 



in Fig. QJi. The state \a) corresponds to an atomic BEC, 
which is coupled to an excited diatomic molecular BEC 
(the state |e)) via a Raman laser pulse Q p ("pump"). A 
single-photon detuning A is the difference between a fre- 
quency of the pump laser and the frequency of transition 
between |o) and |e). The state |e) is, in turn, coupled to 
the ground state of the molecular BEC by a laser pulse 
fid ("dump"). In other words, a pair of atoms from \a) 
is photoassociated by a pump field of molecular Rabi fre- 
quency Q p into a molecule in the excited state |e), which 
is subsequently driven into the ground molecular state. 
This system is intrinsically nonlinear, as it involves merg- 
ing of atoms into diatomic molecules. Its linear counter- 
art is widely used in quantum manipulations: STIRAP 
is well-established technique for transferring popula- 
tion from \a) to \g) in three-level quantum systems. The 
gist of the technique is as following. Intuitively, one may 
think of transferring firstly the population from \a) to |e) 
by means of £l p pulse, and then from |e) to \g). How- 
ever, the state |e) often suffers from spontaneous emis- 
sion. STIRAP counterintuitive sequence of pulses [|[ con- 
sists of overlapping Q Pi d pulses, with ftd going first (see 
Fig. Q~p). STIRAP process in both linear and nonlinear 
A-systems is based on the existence of a dark state which 
is a superposition of only |a) and \g) states (and there- 
fore almost does not suffer from spontaneous emission 
of light, which explains the notation; another notation 
for the dark state is coherent population trapping (CPT) 
state [6]). An instantaneous dark state is determined by 
values of Q p ,d', population transfer happens via the dark 
state as fl Pt d are slowly changed, with |e) state remain- 
ing almost unpopulated. Stability and adiabaticity of the 
system in the dark state is essential for efficiency of pop- 
ulation transfer. In contrast to its linear counterpart, 
the theory of nonlinear STIRAP is not well developed 
yet. This process is very important for BEC physics, as 
it allows to achieve BEC of molecules @ . Mean- field ap- 
proach has shown to be reliable in such systems [|| 0, [l(| • 
Neglecting mean-field collisional interactions and sponta- 
neous emission, the mean-field equations are Q 
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where amplitudes ip a ,e,g ar e normalized as \ip a \ 2 + 
^dV'gl 2 + IV'el 2 ) = 1- The dark state vector is given 
by (up to a phase factor) 
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n d +n c > with = 
8il 2 . Linearization about the dark state 0] gives 



a dynamical system with three eigenfrequencies: ujq — 
0, u± = |[A ± (A 2 + O^) 1 / 2 ]. The frequencies w ,± 
are all real. However, it does not guarantee that the sys- 
tem is always dynamically stable: there might be nonlin- 
ear instabilities. 

We recast the model into the form of a two-degree-of- 
freedom (2 d.o.f.) classical Hamiltonian system, and use 
the classical adiabatic theory and the resonance normal 
forms theory available in [ll|, Eqs. |T]) are equiva- 
lent to Hamiltonian equations of motion of the effective 
classical Hamiltonian 
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Here Xk are canonical momenta, while yk are the coordi- 
nates, being related to the old "variables" (complex num- 
bers ipi) as ip a = xi+iyi, tp e = x 2 +iy 2 , 4> a = x 3 +iy 3 . 
The system has two degrees of freedom only: there exists 
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integral of motion N = xf + yf + 2 (x\ + y 2 + 

Firstly, we consider the case A = (single-photon res- 
onance). Let H = 0. On this manifold, dynamics at 
constant parameters ft p .d is completely integrable. In- 
deed, I 2 = x 2 /y 2 is the additional integral of motion. As 
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X\y\ are 



a result, I 3 = I 2 x 3 + y 3 and I\ 
also integrals of motion. Integrability of the system with 
A = at H = energy manifold is an important finding 
of this paper. In case £l Pt d are changing with time, the 
system may leave the H — energy manifold. However, 
in case initial conditions are such that Ji,3 = 0, then even 
with time-dependent parameters dynamics will be con- 
fined to the initial H = manifold. The particular case 
considered in [|| (where all population is initially in |a)) 
is of this type. Indeed, initially 2:2,3 = 2/2,3 = (there- 
fore, H = 0). As parameters are changed, I\ , I 3 and H all 
remain equal to zero, see Fig. [3J where the example of [6[ 
is presented. The pulses are f2 Pi< j = Sloe~( t ~ tp - d ^ , with 
td < t p . With H = const = 0, we can reduce the system 
to a 1^ d.o.f. classical Hamiltonian. Without loss of gen- 
erality, let I 2 — 0. Therefore, only x\, x 3l y 2 are non-zero 
during the evolution. The equations of motion are x\ = 
fl p xiy 2 , x 3 = ^-y 2: y 2 = -^x 3 ~^x\. Introducing 
new variables <fi and z as x 3 = A= yf\ — x\ cos tfi, y 2 = 

exp[z], the equations of motion 
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correspond now to the classical Hamiltonian 
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FIG. 2: (Color online)Population dynamics and the adiabatic 
parameters A„i, r n i- Tildes over the parameters denote divi- 
sion by 10 (in order to fit into the scale): X = X/1Q. -Pi, 2, 3: 
populations of |a), |e), \g), i.e. x\ + y\ , 2{x\ + y%), 2{x\ + y|) 
correspondingly. Solid lines: Pi, 3; dashed lines: P2 (should 
be zero in the dark state); dot-dashed lines: the adiabatic 
parameters, (a) The same numerical example as in Ref. 0|; 
t p — 3.8, td = 3, Qo = 5. Note the dynamics of P2: at t ~ 3 
it deviates from zero considerably. Parameter r n i suggests 
adiabatic evolution at 1.5 < t < 4.5. Parameter A n i is huge 
(almost two orders of magnitude larger than r n i), suggesting 
nonadiabatic evolution. Note that H, 7i,3 remains zero. (b,c) 
f2o = 100, other parameters are the same as in the previous 
numerical example; (b) Population dynamics. High-frequency 
oscillations are seen in Pi, 2, 3 (c) Dynamics of adiabatic pa- 
rameters A n i , r„i and deviation of Pi from its value in the dark 
state (P?): 8 Pi = Pi - P? . Three stages of nonlinear STI- 
RAP can be defined (I-III), with only the second stage being 
adiabatic. In the first stage (t < 2.5) A n i is large, so system 
acquires some action. At the second stage (2.5 < t < 4.5) A n i 
is small, and population transfer proceeds (with approximate 
conservation of the acquired action) . At the third stage (from 
t « 4.5), adiabaticity is broken again. Predictions of A n i are 
better than that of r n ;. 



where <fi and z are canonically conjugated variables. Fig. 
[TJ; presents phase portraits of this Hamiltonian. There 
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is a stable fixed point which corresponds to the dark 
state: z — | M n+fa ], (j) — 7r - As time increases, 
it goes from z = to — oo. For small-amplitude oscilla- 
tions, linearization around the stable fixed point gives the 
frequency uj — ^(ild^e) 1 ^ 2 ■ Adiabatic evolution requires 
change of this frequency 8u> at one period of unperturbed 
motion to be much less than frequency itself [111, [l3l |: 
\Su\ <C \oj\. Since 5u> « ui * 2ir/u>, we have the following 
criteria: |cj| < lu 2 /2tt, or \Q e tl d + tl e Q d \ < i(f2 d Q e ) 3 / 2 , 
or 



A nl = 2ir\n e n d + n e n d \/{n d n e ) 3 / 2 « 1 



(4) 



It can be seen from Fig. [5]that A n i works sufficiently bet- 
ter than r n i of Q (A n i is typically an order of magnitude 
larger than r n i in this numerical example) . For one-and- 
half d.o.f. Hamiltonian systems, the criterion \Su\ <C \u>\ 
as described above is reliable to analyze adiabaticity. 

In systems with several degrees of freedom, situation is 
much more complicated. Even stability at fixed parame- 
ters is highly nontrivial issue. Indeed, let us now consider 
the case A^O (or A = 0, but H ^ 0), where we cannot 
utilize the trick with the energy manifold, so the system 
with fixed parameters Sl P:d has 2 d.o.f. We make several 
canonical transformations to reduce the Hamiltonian sys- 
tem to a 2 d.o.f. one, and shift the o rig in to the point 
corresponding to the dark state (see 14 1). We get the 
Hamiltonian in new variables gi,2,Pi,2 
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H = n p [q 1 +p 1 +q 2 +p 2 ]p 1 — —qiq 2 — —piP2~^ (<7i+Pi)- 

One may transform the quadratic part of the Hamilto- 
nian to its normal form: sum of two linear oscillators, and 
neglect the cubic part, obtaining therefore two uncoupled 
oscillators (with frequencies lo±). However, such straight- 
forward approach is dangerous at low-order resonances 
between tu±. The 1:1 resonance happens at A = 0, while 
1:2 resonance at A = y/il d il e /8. Consider firstly 1:1 
resonance (A = 0). We introduce a parameter 9 as 



Q d = fl e cos9, il p = £l e sin#/2v2, 



and divide the Hamiltonian ((5 



by f2 e 



1 > cos 9 > 0, 



H = s ^(q 2 



ji - ^-qiq 2 ~ p>\P2- 
We transform the quadratic part of the Hamiltonian 
to the normal form by means of linear transformations 
and obtain (retaining old notations for new variables) 

1 



H = H 2 + H 3 = - V cos 0(p 2 ?i - Pi q%) + H 3 , (6) 

where H 3 consists of cubic terms (for 1:1 resonance, it is 
convenient to use series in p 2 qi — Piq 2 , Pi + p\, q\ + <?!■ 
Hamiltonians which depends only on such combinations 
of variables are also called normal forms). We need to get 
rid of the cubic terms in the Hamiltonian. To this end, we 
fulfill a nonlinear canonical transformation using a gen- 
erating function K 3 which is a homogenous polynomial 
of the third order in old coordinates and new momenta: 
K 3 = aq\q 2 + f3P 2 q 2 + zP 2 q 2 + lq \ + fPiP 2 qi- For the 
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FIG. 3: (Color online) Instability in 1:1 resonance. Dynamics 
of population of the excited state P2 (which should be zero 
in the dark state) is shown. Parameters: Q p — 2.0, fid = 1-0 
(cos# ~ 0.17 < 0.2, Q p > \/3^d)- Instability develops on 
long timescales: we start with small deviation from the dark 
state, but the deviation grows and reach very large values. 
With cos# > 0.2 (not shown), initial deviation do not grow 
much: it slowly oscillates about the equilibrium of the effec- 
tive potential of (7). The curve P2{t) is filled by oscillations 
with frequency ^V^dSle which are not seen in this scale. In- 
set shows two effective potentials of (0 at J = 0.001 and 
two different values of 9; solid line with cos# = 0.17 < 0.2 
corresponds to the unbounded motion (instability). 



coefficients a,/3,7, /, z in K 3 , we get a system of linear 
equations; the transformation defined by K 3 kills all cu- 
bic terms, while quadratic part remains the same; at the 
same time, quartic terms emerge: we get the Hamilto- 
nian H = iVcos6»(P 2 (3i - P1Q2) + H 4 + .., where H 4 
contains the quartic terms. Generally, it is not possi- 
ble to kill all quartic terms [ll[. We need to transform 
i?4 to the normal form. However, we can proceed in a 
simpler (but equivalent) way. We change to polar coor- 
dinates in Qi,Q 2 plane making a transformation Qi — 
i?cos£, Q 2 = i?sin£, Pi = Pcos£- P 2 = 

Psin£ + J s , and then average over the (fast) angle £. 
We obtain the following effective Hamiltonian: 



1 Jul 
p — _ j <L 

r s 



A, 



-(P 2 



R 2 



) + R 2 A R 



(7) 



where A P = 10cos2 9 + 3 1 cose + 4 i A R = 5cos<9 - 1, k 6 = 
, 81 " e In this Hamiltonian, J is constant. Therefore, 

\/2Vcos e ' ' 

it is an integrable system for a pair of canonically conju- 
gated variables P, R. It is not difficult to understand its 
dynamics: ([7]) is merely a Hamiltonian of a particle in the 
potential U(R) = J 2 ^f + ArR 2 (indeed, only the expres- 
sion in the square brackets in ((7|) is important for dynam- 
ics). The effective potentials U(R) are shown in the inset 
of Fig. ([3]). For cos 9 > -g, the system (J7]) has a single fixed 
point, while for cos9 < | there are no fixed points. In 
the latter case (corresponding to S~l p /Q d > a phase 
point initially placed close to R = will slowly move 
from the origin to large values of R. The dark state cor- 
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responds to R = 0, therefore this case implies dynamical 
instability due to 1:1 resonance. The instability is slow 
(see Fig|3]), nevertheless it is physically important: it de- 
velops on the same timescales as that used in [7|. We 
checked that crossing the critical cos 8 — ^ value results 
in onset of dynamical instability in accordance with the 
effective potential ([7]). Physically, cos 9 < -| means that 
most of the population is in \g) state. In this region, the 
dark state is unstable: deviations from it slowly grows 
with time, populating |e) state considerably which would 
lead to spontaneous emission in real experiments. To be 
more precise, deviations from the dark state firstly os- 
cillate in time with frequency ^y/Tt^Tl and small ampli- 
tude, and then the amplitude of oscillations slowly grows. 
Thus, one may see periodic bursts of spontaneous emis- 
sion each time the deviations reach their maxima (i.e., 
with a time period T ~ l/vTW^e)- 

Let us now turn to 1:2 resonance. Usually, its analysis 
is much simpler than that of 1:1 and 1:3 resonances, as 
it involves only cubic terms. Unfortunately, in our case 
this resonance is also degenerate. In 1:2 resonance, we 
can bring the Hamiltonian to the resonance normal form 

H res = wifli + uj 2 R 2 + kR{ /2 R 2 sin^i + 2<j> 2 ), (8) 

where (Ri,4>i) are symplectic polar coordinates (action- 
angle variables of H 2 : P% = V^Ri sin fa, Qi — 
V2i?i cos fa), k is some coefficient, the frequencies u>i and 
lu 2 are in 1:2 resonance: lui + 2lu 2 = (so, W1W2 < 0). The 
resonance normal form is obtained by averaging out all 
terms in H 3 except the resonance one (which depends on 
the resonance phase 7 = fa + 2fa) pT| . The resonance 
normal form is integrable. Provided certain conditions 



of generality are fulfilled, the equilibrium of the original 
system (i.e., the dark state) is stable or unstable simul- 
taneously with the equilibrium of the resonance normal 
form (i.e., R\^ 2 = 0). For 1:2 resonance, a condition of 
generality is k ^ 0' with uj\lo 2 <C and k 0, the equilib- 
rium is unstable [1_J. In our case, n = 0, therefore this 
source of instability is absent in the system. Still, the 
degeneracy k — does not guarantee stability but des- 
ignates that possible instabilities, if exist, are very slow. 

It is important to emphasize that stability and adia- 
baticity are two different issues. Adiabatic invariance in 
systems with several degrees of freedom is also remark- 
ably nontrivial [!_}. Main difficulties come from passage 
through resonances; in case our system stays away from 
low-order resonances as parameters are changed, we can 
consider it as two decoupled oscillators and generalize the 
adiabatic condition @ straightforwardly see [l4T |. 

Here, we analyzed low-order resonances of the non- 
linear A system; found nonlinear instabilities of the dark 
state; found cases of complete integrability and improved 
adiabatic conditions for nonlinear STIRAP. The sug- 
gested method is generalized to systems with mean-field 
collisional interactions in [14j |. Briefly, theory of nonlin- 
ear STIRAP is developed. 
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